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One possible fate of an accreting white dwarf is explosion in a type la supernova. 
However, the route to the thermonuclear runaway has always been uncertain owing to the 
lack of a convective model consistent with the Urea process. 

We derive a formalism for convective motions involving two radial flows. This formalism 
provides a framework for convective models that guarantees self-consistency for chemistry 
and energy budget, allows time- dependence and describes the interaction of convective 
motions with the global contraction or expansion of the star. In the one-stream limit, we 
reproduce several already existing convective models and allow them to treat chemistry. 
We also suggest as a model easy to implement in a stellar evolution code. 

We apply this formalism to convective Urea cores in Chandrasekhar mass white dwarfs. 
We stress that in degenerate matter, nuclear reactions that change the number of electrons 
strongly influence the convective velocities. We point out the sensitivity of the energy 
budget on the mixing. We illustrate our model by computing stationary convective cores 
with Urea nuclei. We show that even a very small mass fraction of Urea nuclei (10 -8 ) 
strongly influences the convective velocities. 

Finally, we present preliminary computations of the late evolution of a close to Chan- 
drasekhar mass C+0 white dwarf including the convective Urea process. 



1. Prelude to thermonuclear explosions 

When a CO white dwarf accretes matter, its centre gets denser and hotter. At some 
point, the C-C burning ignites mildly in the centre. When the radiative luminosity can 
no longer get rid of the heat produced, a convective core forms and grows as the burn- 
ing releases more and more energy. During the burning process, Urea pairs such as 
23 Na- 23 Ne are produced. These release neutrinos through emission and capture of elec- 
trons preferentially around Urea shells. The net amount of energy released and the change 
in the electron fraction at the time of the explosion have always been uncertain due to 
the lack of a convective model self-consistent with the energy and chemistry budgets [ 
EI21IIIlSISlini0IEllI3[IIIl- Former theoretical attempts [ El Q2 [13] to describe the 
interplay between convection and chemistry have all left open the question of energy con- 
servation. 



A better understanding of the phase immediately preceding the thermal runaway pro- 
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vides an essential link between the core evolution in the progenitor and the supernova 
explosion. Knowledge of the physical and chemical profile at the time of the explosion 
may also provide new insight in the physics of the explosion itself. Indeed, the neutron 
excess at the time of ignition is a crucial parameter of the explosion, and it is strongly 
affected by the convective Urea process. 

2. A model for the convective Urea process 

Except for the work of Iben [ |1J E], the convective Urea process has been so far 
neglected in presupernova computations. This is because we lack a proper energy equation 
in a convective region. Indeed, far from the Urea shell, out of equilibrium nuclear heating 
overwhelms neutrino losses [ E] . Hence the net heating crucially depends on the actual 
chemical mixing. The amount of convective work necessary to carry the electrons is 
uncertain [ El E] Last but not least, in degenerate matter chemistry controls the 
buoyancy and thence the convective velocities. Fortunately, the problem is well known 
and solved in the radiative state. Our basic idea is to mimic convection by considering 
the interaction between two radiative columns of fluid side by side [I14j. 

The resulting convective model guarantees energy conservation. Besides, it is 
time-dependent, non-local, allows for different rising and descending convective velocities, 
describes horizontal chemical inhomogeneities and treats the interaction of convection 
with the mean flow. We have checked carefully different limits of our model against 
corresponding existing convective models and found good agreement. 

Although the full model would be rather difficult to implement in a convective core, the 
one-stream, local, time independent limit yields a straightforward extension of the mixing 
length theory fully consistent with chemistry and energy budget. The main difference with 
the mixing length is the dependence of the convective velocity on the chemical gradients. 

3. Application to degenerate CO convective cores 

We use our convective model to compute convective cores with a homogeneous total 
abundance for the Urea pair A=23. Indeed this pair operates at mass density 1.7 10 9 g cm - " 
and is hence one of the main actors at the onset of C burning which occurs around density 
2 10 9 g cm -3 and temperature 3 10 8 K. 

3.1. Stationary models 

We first used the full two-stream model to compute stationary convective cores. We 
varied the mass fraction of Urea matter Xy from 10 -12 to 10 -3 . Strong modifications of 
the convective profile occur for mass fractions higher than 10 -8 as Fig. Qshows. Analytical 
arguments on the existence of a stationary convective core confirm this low threshold. 

The convective velocities are strongly reduced by the presence of Urea matter. Indeed, 
the velocity of a convective blob is proportional to its difference in density with the 
surrounding medium. Besides, the density of degenerate matter depends mainly on the 
electron number density. This means that convection in degenerate matter crucially 
depends on electron emissions and captures. These simulations also show that the 
one-stream, local limit is a good approximation to the two-stream model in this context. 
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Figure 1. Variation of the stationary convective profile for different mass fractions X\j of 
the Urea pair A=23. 

3.2. Stellar evolution models 

We have started to implement a one-stream, stationary limit of the two-stream model 
in P. Eggleton's stellar evolution code [ llfij . This has already required significant changes 
to the code to handle chemical gradients with sufficient accuracy. 

First, it appears that the backreaction of the chemical gradients on the convective 
velocity, and hence the mixing, causes serious convergence problems. We momentarily 
overcome these problems by artificially lowering the mixing at the outer edge of the 
convective core. Simulations with X\j = 10~ 7 confirm the reduction of the convective 
speed which now takes place only upward of the Urea shell. 

Second, for higher Urea mass fractions, our code breaks down when the outer edge of 
the convective core reaches the Urea shell. At this point, the convective state changes on 
short time scales due to the gradient of the Urea composition. 

4. Conclusion and future prospects 

We have developed a new convective model suitable for studying the convective Urea 
process. A first application to stationary convective cores suggests that electron emissions 
and captures control convection in degenerate matter. 

We hence expect to uncover new results on the initial conditions for the explosion of 
type la (and possibly core-collapse supernovae). In particular, we shall be able to predict 
the degree of neutronisation at the time of the thermal run-away, a crucial parameter for 
explosive nucleosynthesis 

However, we still need to implement a fully time- dependent version of our convective 
model in a stellar evolution code including an extended chemical network. 
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Finally, our model may well prove to be relevant for other fields of stellar evolution 
such as the He flash, the fate of O, Ne and Mg cores and AGB stars. 
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